Modified Dispersion Relations lead to a finite Zero Point Gravitational Energy 
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We compute the Zero Point Energy in a spherically symmetric background distorted at high 
energy as predicted by Gravity 's Rainbow. In this context we setup a Sturm-Liouville problem with 
the cosmological constant considered as the associated eigenvalue. The eigenvalue equation is a 
reformulation of the Wheeler- De Witt equation. With the help of a canonical decomposition, we 
find that the relevant contribution to one loop is given by the graviton quantum fluctuations around 
' the given background. By means of a variational approach based on gaussian trial functionals, we 

find that the ordinary divergences can here be handled by an appropriate choice of the rainbow's 
• functions, in contrast to what happens in other conventional approaches. A final discussion on the 

' connection of our result with the observed cosmological constant is also reported. 

oo 

I. INTRODUCTION 

u\ 

'■ The idea of promoting General Relativity to a quantum level, together with a solution of the Cosmological Constant 
5-H Problem, is one of the big challenges of our century. Indeed a satisfying Quantum Gravity theory does not exist yet and 
I the enormous gap of 10^^" orders of magnitude between the predicted theoretical value of the cosmological constant 

and the observed one has not yet found a compelling explanation. If the quantum description of Nature is appropriate 
I ] for every force, it should be applicable even to the gravitational force described by General Relativity. But perhaps, 
K>" ' General Relativity as it stands, requires a change In this respect, various proposals on how the fundamental aspects of 
special relativity can be modified at very high energies have been done. Among these proposals, particularly promising 
i appears to be the one known as Doubly Special Relativity (DSR) 1]. One of the characterizing DSR effects is that the 
' usual dispersion relation of a massive particle of mass m is modified into the following expression 

g '■ E^gl [E/Ep) - p^gl (E/Ep) = m^, (1) 

where gi (E/Ep) and g2 {E/Ep) are two arbitrary functions which have the following property 



X 



lim gi (E/Ep) = 1 and lim (E/Ep) = 1. (2) 

E/Ep^O E/Ep^a 



^ . Thus, the usual dispersion relation is recovered at low energies. Of course the first ideas of DSR were minted for fiat 
^ space. However, nothing forbids to consider a curved background and therefore to enter into the realm of General 
Relativity. From this point of view, Magueijo and Smolin[2] proposed that the energy-momentum tensor and the 
Einstein's Field Equations were modified with the introduction of a one parameter family of equations 

G^, (E) - 8^G (E) T^, {E) + g^,K (E) , (3) 

where G (E) is an energy dependent Newton's constant, defined so that G (0) is the low-energy Newton's constant. 
Similarly we have an energy dependent cosmological constant A (E) leading to the rainbow version of the Schwarzschild 
line element 

^2 A 2MG(0)\ dP P .^^2 ■ 2n^j.2\ 

ds^ = - [1 — — ■ + -, r h — ; r [dO^ + sm^ dd(^^\ . (4) 

V r ) a\[E/Ep) U-^^^f^Agl^E/Ep) gliE/Ep)"^ ^> 



* Electronic address: Remo. GarattiniOunibg. it 

t Electronic address: Gianluca.Mandanici@unibg.it, 



2 



Since the functions gi (E/Ep) and (72 {E/Ep) come into play when the energy E is comparable with Ep, it is likely 
that they modify the UV behavior in the same way as Generalized Uncertainty Principle and Noncommutative 
Geometry (NCG) do, respectively. If the effect of Generalized Uncertainty Principle and NCG is to modify the 
Liouville measure cPxcPk from one side, and to introduce a granularity from the other side, the rainbow metric 
should be able to introduce a natural UV regulator hidden into the arbitrary functions gi {E/Ep) and 52 [E/Ep). 
An encouraging partial answer has been obtained in Ref. [3] , where an application of Rainbow's Gravity to black hole 
entropy computation has been considered. In that paper the UV regulator, namely the brick wall, has been eliminated 
with the help of the following choice of gi {E/Ep) and 52 {E/Ep) 

^^(^/^-^-expf-AV (5) 



32 {E/Ep) Ep 

An interesting test to see Rainbow's Gravity at work again, should be the computation of Zero Point Energy (ZPE). 
Nevertheless, we have to remark that any computation of ZPE leads to a regularization and subsequently to a 
renormalization process in order to have finite physical quantities. Therefore the purpose of the paper is to show 
that with the introduction of appropriate Rainbow's Gravity functions it is possible to overpass the renormalization 
problem. Of course this proposal does not represent a complete cure to have a finite theory of Quantum Gravity but 
rather it suggests how the modification of some basic principles like the introduction of an energy dependent metric 
can lead to unexpected results such as the avoidance of a renormalization scheme. In ordinary gravity the computation 
of ZPE for quantum fluctuations of the pure gravitational field can be extracted by rewriting the Wheeler-DeWitt 
equation (WDW)[i] in a form which looks like an expectation value computation[5]. We remind the reader that 
the WDW equation is the quantum version of the classical constraint which guarantees the invariance under time 
reparametrization. Its original form with the cosmological term included is described by 



{2>i)G,,kin'^TT''' -^{^R-2A) 



* = 0. (6) 



If we multiply Eq.([5]) by ^* [gij] and functionally integrate over the three spatial metric gij, we can write^Q 



1 JV[g,,]^* [.g.,]/srf3xAs*[.9y] 



4- 



A 



V JV[g,,]^*[g,,]^[g,,] V {^\^) k' 

where we have also integrated over the hypersurface S and we have defined 



(7) 



y = / d^xv^ (8) 



as the volume of the hypersurface S with 

As = (2k) G^jkin^'n'^' - ^'R/ (2k) . (9) 

In this form, Eq.Q can be used to compute ZPE provided that A/k be considered as an eigenvalue of As, namely the 
WDW equation is transformed into an expectation value computation. In Eq.([B]), Gijki is the super- metric, vr*^ is the 
super-momentum, ■^i? is the scalar curvature in three dimensions and A is the cosmological constant, while k = SvrG 
with G the Newton's constant. Nevertheless, solving Eq.([7]) is a quite impossible task, therefore we are oriented to 
use a variational approach with trial wave functionals. The related boundary conditions are dictated by the choice of 
the trial wave functionals which, in our case, are of the Gaussian type. Different types of wave functionals correspond 
to different boundary conditions. The choice of a Gaussian wave functional is justified by the fact that ZPE should 
be described by a good candidate of the "wacMwm state'' . To fix the ideas, a variant of the line element Q will be 
considered 

ds^ = -N^ {r) + + {d9'' + sin2 edcf>^) , (10) 

(l-^JL)^g^^{E) 9l{E)^ 

where N is the lapse function and h (r) is subject to the only condition b{rt) = rt- Metric (1101) . will be our cornerstone 
of the whole paper which is organized as follows. In section|lll we derive the Hamiltonian constraint in presence of the 



^ See also Ref. 0| for an application of the method to a / (_R) theory. 
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background (fTO)) . in section Hill we compute the ZPE of quantum fluctuations around the background ^TU\\ and with 
the help of an appropriate choice of the functions gi {E/Ep) and 52 (E/Ep), we will show that the UV divergences 
of ZPE disappear. We summarize and conclude in section ITVl Units in which h ^ c = k = 1 are used throughout the 
paper. 



II. THE HAMILTONIAN CONSTRAINT IN RAINBOW'S GRAVITY 



In order to use Eq.^ for the metric ^TU\\ . we need to understand how the WDW modifies when the functions 
gi (E/Epi) and 52 (E/Epi) distort the background. It is therefore necessary to understand how some basic ingredients 
change under the transformation of the line element ()10|) . The form of the background is such that the shift function 



N' = -Nu' = g^' 

vanishes, while N is the previously defined lapse function. Thus the definition of Kij implies 



9ij 
2N 



(11) 



(12) 



where the dot denotes differentiation with respect to the time t and the tilde indicates the quantity computed in 
absence of rainbow's functions g\ (E) and g2 (E). For simplicity, we have set Ep — 1 in 171 (E/Ep) and g2 (E/Ep) 
throughout the paragraph. The trace of the extrinsic curvature, therefore becomes 



K = g'^K,,^g, [E) K 
and the momentum tt*-' conjugate to the three-metric gij of E is 

2n^^ > 52 {E) 



(13) 



(14) 



Now, we have enough information to define the WDW equation for a background described by ([TU]). From Eq.® we 
find that = becomes 



V9 



2Kg2 [E) 



2Ac 



* 0, 



where we have used the following property on R 

R^g'^R^^^gl {E) R 

and where 

1 



G 



ijkl 



{gikgji +giigjk - gijgu) 



G 



ijkl 



2^ w..... . ...... ....... g^^j^y 

Therefore, in presence of Rainbow's Gravity, we find that Eq. ^ becomes 



(15) 



(16) 



(17) 



gl {E) (* k ^'^As 



where 



Ay 



V 



{2n\ 



Ac 

K 



.9? (E)^ 



gUE) 



VgR 



(2«)<?2 (E) 



(18) 



(19) 



We can gain more information if we consider gij = gtj + hij, where gij is the background metric and hij is a 
quantum fluctuation around the background. Thus Eq. ljlSp can be expanded in terms of hij. Since the kinetic part of 
As is quadratic in the momenta, we only need to expand the three-scalar curvature / (Px^/g^R up to the quadratic 
order. However, to proceed with the computation, we also need an orthogonal decomposition on the tangent space of 
3-metric deformations 0, Q: 
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(20) 



The operator L maps into symmetric tracefree tensors 



(21) 



hjj is the traceless-transverse component of the perturbation (TT), namely 



V'h. 







(22) 



and h is the trace of hij. It is immediate to recognize that the trace element u = /i — 2 (V • ^) is gauge invariant. It is 
straightforward to see that the gauge invariant decomposition (|20p does not change, when we consider the rainbow's 
metric ([TU]). Therefore, following the results of Ref.Q, we can use the final expression 



(2) 



1]/" 



(2) 



V 



V 



K 



(23) 



Note that in the expansion of J-^ d^x^fgR to second order in terms of hij, a coupling term between the TT component 
and the scalar one remains. However, the Gaussian integration does not allow such a mixing which has to be introduced 
with an appropriate wave functional. Extracting the TT tensor contribution from Eq. (|18p . we find 



9l{E) 



d-'xy gG 



-lijkl 



{2k) 



9UE) 



{x,x) 



ijkl 



{2k) 92 {E) 



A^K^ {x,x) 



ijkl 



(24) 



The origin of the operator comes from 



(25) 



which is the modified Lichnerowicz operator where A^is the Lichnerowicz operator defined by 

{AlH) = Ah,j - 2R,kjih''' + R,khh + R,kh\ A = -VVa- 



(26) 



Qijki represents the inverse DeWitt metric without the y/g factor and all indices run from one to three. Note that 
the term 



(27) 



disappears in four dimensions when we use a background which is a solution of the Einstein's field equations without 
matter contribution. The symbol in Eq. (p4)) means that we have rescaled every piece in Eg.p^. evaluated at 
second order. Moreover, although the expression of As explicitly shows how globally changes the operator As, when 
we consider the following eigenvalue equation 



(Ar.-), 



E^h. 



(28) 



we find that 



E' 



9i{E) 



(29) 



in order to reestablish the correct way of transformation of the perturbation. Then, the propagator (x, x)^^i.i can 
be represented as 



iakl 



E 



^l:^"(^)^i^^"(^) 

2A(r)5|(i?) ' 



(30) 
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where C^^) the eigenfunctions of . r denotes a complete set of indices and A (r) are a set of variational 

parameters to be determined by the minimization of Eq. (j24[) . The expectation value of is easily obtained by 
inserting the form of the propagator into Eq. (I24p and minimizing with respect to the variational function A (t) . Thus 
the total one loop energy density for TT tensors becomes 



A 
8^ 



(31) 



The above expression makes sense only for Ef (r) > 0, where E.i are the eigenvalues of . With the help of Regge 
and Wheeler representation 10], the eigenvalue equation can be reduced to 



rf2 ^(/ + 1) ■ 
— + ^2 + ('^) 



dx 



E? 



9l{E) 



1,2 , 



(32) 



where we have used reduced fields of the form fi (x) — Fi (x) /r and where we have defined two r-dependent effective 
masses m\ (r) and m| (r) 



nl{r)^^(l-'-^)+^b'{r)-^b{r) 

-iW = ^(i-^ 

In order to use the WKB approximation, from Eq 



(rEEr(x)). (33) 
we can extract two r-dependent radial wave numbers 



Elnl 1(1 + 1) 



gUE) 



(34) 



III. ONE LOOP ENERGY IN AN ORDINARY SPHERICALLY SYMMETRIC BACKGROUND 



It is now possible to explicitly evaluate Ea.([3T|) in terms of the effective mass. To further proceed we use the W.K.B. 
method used by 't Hooft in the brick wall problem [lll| and we count the number of modes with frequency less than 
uji, i = 1,2. This is given approximately by 



'(E,)^ / ,^,{l,E,)i2l + l)dl, 



where Vi (/, Ei), i = 1, 2 is the number of nodes in the mode with (/, Ei), such that {r = r {x)) 

u^{l,E,) = - / dx^k1[r,l,Ei). 



(35) 



(36) 



Here it is understood that the integration with respect to x and /max is taken over those values which satisfy 
kf (r, Z, Ei) > 0, i = 1, 2. With the help of Eqs. p5|36p . Eq.dSI]) leads to 



A 
8^ 



1 2 r.+ca 

-E/ E,g^{E)g2{E) 

^=l 



d-g {E, 
dE, 



■dE, 



(37) 



This is the graviton contribution to the induced cosmological constant to one loop. The explicit evaluation of the 
density of states yields 



dEi 



dE, 



+ 00 pi 

dx 



^kHr,l,E)dE, \gl{E) 



E 



2 /■+°° , 2^ _ 
7ri_ ^"^^ dE, \gl[E) 



mf (r) 



E} A f+°^ , , d f Ef 



9l (E) 



-mfir) = — dxr'—[-^-mf{r) 



dE, \gl{E) 



(38) 
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Plugging expression (|38l) into Eq. (|37)) and dividing for a volume factor, we obtain 

S^ = -^T.l ^^9. iE) (E) - ™f (O) dE., (39) 

where E* is the value which annihilates the argument of the root. In the previous equation, we have included an 
additional 47r factor coming from the angular integration and we have assumed that the effective mass does not 
depend on the energy E. To further proceed, we can see what happens to the expression (j39p for some specific forms 
of gi {E/Ep) and 32 (E/Ep). One popular choice is given by 

gi{E/Ep) = l-7]{E/Ep)'' and g2iE/Ep) = l, (40) 

where 77 is a dimensionless parameter and n is an integer[12]. Nevertheless, the above choice does not allow a finite 
result in Eq. (p9| and therefore will be discarded. Thus the choice of the possible forms of gi [E/Ep) and 32 {E/Ep) 
is strongly restricted by convergence criteria. We have hitherto used a generic form of the background. We now fix 
the attention on some backgrounds which have the following property 

TOq {^) = "^2 — —iril (r) , Vr G {rt,ri) . (41) 



For example, the Schwarzschild background represented by the choice b (r) = rt = 2MG satisfies the property (|4T|) 
in the range r g [rt,5r(/2]. Similar backgrounds are the Schwarzschild-de Sitter and Schwarzschild- Anti de Sitter. 
On the other hand other backgrounds, like de Sitter, Anti-de Sitter and Minkowski have the property 

™o ~ '^2 ('") = "^1 (^) 1 G ('"t; • (42) 

In case condition (|4T|) holds, Ea. (p9)) becomes 

A _ J_ 

8^ ^ ^3^ 

where 



(/++/-), (43) 



and 



/■°° d f E"^ \^ 

1+ = (Egi {E/Ep) 52 [E/Ep)) — [-J^^J^ + < (0 j dE (44) 

/•°o d / E'^ \^ 

I = I (Eg, iE/Ep) 52 (E/Ep)) - (^-^^^ - (r) j dE. (45) 



Instead, in case condition (|42l) holds, Eq. ([M|) becomes 

A 2 



BttG 37r 



(46) 



We begin to look at Eq. p5|) . It is immediate to see that integrals /+ and /_ can be easily solved for a very particular 
choice. Indeed, if we set 

g^HE/Ep)^g,{E/Ep) (47) 

we find that and /_ take the form: 



E V d f E \ // E 



and 



E V d f E \ If E 
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The above integrals can be easily evaluated using the auxiliary variable 

E/Ep 



z (E/Ep) = 



92 {E/Ep) 



so that Eq. (|39|) becomes: 



A 
8^ 



E% 



' \/ — x^dz + J z^ \/ z^ + x^dz \ , 



(50) 



(51) 



where z^o = Iure^oo z{E/Ep) and x = \Jrn^ (r) /-Ep. The integrals involved in Ea. ([ST|) can be calculated straight- 
forwardly being: 



and 



/i,^(z) = y zV^^ - x^dz = g {z {2z^ - x^) Vz^ -x^-x^ log 2 (z + ^/ z^ - | 
/2,x(z)= y" zV^^ + x'^dz ^^{z {2z^ + x^) Vz^ +x^-x^ log 2(z + VIM~J?2^ I 



Thus we get the final expression: 

A 



Ej, 



(52) 



(53) 



(54) 



for the case with x < Zqo, and the expression: 

A 



E'p 



SttG 



{h.xiZoo) - h.xiO)} , 



(55) 



for the case with z^o < x. In particular for the class of rainbow functions that satisfy the condition Zoo = we get a 
vanishing cosmological constant: 



A 
8^ 



0. 



(56) 



Although very appealing, the result (1551) presents the unpleasant feature of being always negative, even in the region of 
space where we would expect a positive cosmological constant. Moreover it is independent on the choice of g{E/Ep) 
provided that this last one can guarantee the convergence of the integral and the absence of imaginary factors. For 
these reasons we are led to investigate other forms of gi{E/Ep) and g2{E/Ep) even if they have less symmetry 
with respect to proposal (|47l) . The only restrictions we have are the low energy limit ^ and the the convergence 
requirement for the integrals p4)) and (l45t . To do calculations in practice a useful choice is the following 



" E"^ E'^ 
gi [E/Ep) = ^^i-_exp(-a-2-), 



92 {E/Ep) 



a > 0, /3i e 



(57) 



The use of a "Gaussian" form is dictated by the possibility of doing a comparison with NCG models. Indeed, in 
Ref. Is*! the authors have considered a distortion induced by an underlying NCG on the counting of states. Basically, 
one finds that the number of states is modified in the following way 



d^xd^k d^xd^k 
dn = 5- =^ dri j = 5- exp 

(2-)' {2^f 



2 



M) 



1,2, 



(58) 



where the UV cut off is triggered only by higher momenta modes > l/Vd which propagate over the background 
geometry. Then the induced cosmological constant becomes 



A 



1 



SttG 611^ 



+00 



(cj2 _ ml (r))^e-^ 



+ 00 



^J{Lu^ + ml (r))^e-^("'+"o(''')da; 



(59) 
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It is immediate to see the analogy with the choice ((57|) . However Eq. (|59)) leads directly to a positive induced cosmo- 
logical constant, while Eq.(|33]) needs an appropriate choice of gi (E/Ep) and 52 (E/Ep) to induce a positive part. 
After choice (f57| , the graviton contribution terms (|44|) and (|45|) become 



and 




In Appendix (fS| . we explicitly compute the integrals (j60p . (|6ip for every n. In order to motivate choice (j57p . we have 
to observe that the case with n = leads to a negative value of A/SttG for every kind of background as one can see 
from Eq. (|43|) . Thus, it is necessary a correction on the pure Gaussian choice in such a way to have a possible change 
of sign in A/SvrG. For our purposes, it is sufficient to discuss the case with n = 1 and n = 3. We begin with n — 1. 



A. Example a) 



After integration, for n = 1, Eq.(j43|) can be rearranged in the following way 

A A , „ . 1 



inGEj, 



— cosh — — 
a \ 2 



ax 



-B — \- smh [ax ) H — cosh (ax ) + 



x^- 



— ) e"^' erf (^^^a;) 
2a ' 



(62) 



where, again, x ~ \J mp (r) jE'p , /3i = /? and where Kq (x) is the Bessel function and erf (x) is the error function. It 
is clear that for every choice of the couple (a, /3) there exists a curve with a different behavior. Therefore to fix ideas, 
we will fix the Gaussian factor a to the same one proposed by the NCG setting of Eg. ((5^1 . Before doing this, it is 
useful to compute the series expansion for small and large x. For large x one gets 



A 



(2/3^3/2 ^ ^^2) ^ g^^5/2 ^ 3^^3 



SnGEj, 



167r2aii/2a 



1287r2 



Q,15/22;3 



while for small x we obtain 



A 



SnGEj, 



It is straightforward to see that if we set 



47r2a9/2 



0(a 



(63) 



(64) 



(65) 



then the linear divergent term of the asymptotic expansion (1551) disappears and Eq. (p^ vanishes for large x. Plugging 
Ea.(|65|) into expansion ([64)) . we obtain 



SirGEj, 



(66) 



which means that for x = 0, the induced cosmological constant never vanishes and therefore cannot be a good candi- 
date to reproduce the Minkowskian limit. Indeed, we have to recall that the variable x expresses the curvature of the 
background through the shape function b (r), which for Minkowski vanishes. On the other hand, the vanishing of ex- 
pression (|63)) and consequently Eg. ((62)) for x — 00 offers a good candidate for large distances estimates. Alternatively, 
by imposing that 



(67) 
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the expression ([64)1 vanishes for small x, while for large x, the leading term becomes 

A (Stt - 8) a; 



iirGE' 



(68) 



This means that Ea.(|62p diverges towards negative values. It is straightforward to see that we cannot simultaneously 
fix both the conditions ([55)) and (1571) for the same a in order to have a vanishing expectation value of A/SttG for small 
and large x unless we consider different values for a for the different behaviors. The idea is to find a point where a 
transition from one parametrization to the other one exists. To begin we have to observe that if we fix one couple of 
parameters to 

= (^ = A\f^' (69) 

4 6 \ TT 

where ai has the same value of the numerical factor appearing in Ea. (|5n| and the second couple with generic values, 
one discovers multiple roots where a smooth transition from one parametrization to the other one can happen. This 
is illustrated in FiglUwhere the couple (|M| together with some generic values of the couple satisfying condition ((55)) 




FIG. 1: Plot of A/SttG as a function of the scale invariant x. Choosing parametrization (|69p . we obtain the vanishing of A/SttG 
when a; — >■ 0. The other curves satisfy condition J65]) for different values of a, with a ^ ai. It is visible the presence of multiple 
roots. 



are shown. It is visible the presence of multiple roots. It is also immediate to see that there exists one and only one 
transition point which can be found by imposing the existence of a tangent point between the curves parametrized 
by the values in (|5^ and the curves parametrized by 

a2, 13 = -^, (70) 
where 0^2 is to be determined. We end up with the following choice 

= .7744164292, l3 = ^^^^y-, (71) 

where the common point is located in x = 1.818231873 as shown in Figl2lThis choice corresponds to the following 
setting 



51 (E/Ep) = cxp(-ai 1^) (l - 7^3^) , 92 (E/Ep) = 1 0<x<l 



818231873 



gi{E/Ep)^exp{-a2^)(l-^^§^) , g2{E/Ep)^l x > 1.818231873 



(72) 
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FIG. 2: Plot of A/SttG as a function of the scale invariant x. For qi = 1/4, there exists a2 — .7744164292 where a smooth 
transition between the two asymptotic behaviors is possible. The transition appears for x = 1.818231873. 



The setting (|72t allows the expression p3)) to have finite values for every kind of background of the spherically 
symmetric type. Let us apply our result to the Schwarzschild background. In terms of the variable a;, we find that 



1^ r > 2MG 




2MG ■ ^''^^ 



Its behavior is 



oo 



when M^O for r ^ 2MG 



^ ~^ I when M ^0 for r> 2MG ' ^'^^^ 
while 

when M oo for r^2MG , . 

oo when M — oo for r > 2MG ' ^ ' 

The situation with Af — > oo describes a wormhole incorporating the whole universe which is not a physical situation, 
while for M — >■ we approach the Minkowski limit which should predict a vanishing induced cosmological constant. 
Note that for both setting ([55)1 and (|67p , we find that the whole behavior can be summarized by the following double 
limit 

A (r) A (r) 

lim lim — — ^ hm lim (76) 

M^ar^2MG SttG r—>2MG M^O SttG 

suggesting that a sort of non-commutativity emerges in proximity of the throat. Therefore, when we adopt the 
parametrization (|72p . the Minkowskian limit is recovered for every value of M. Turning now to the case of Ea.(|46p. 
we find that it is possible to have only one parametrization to obtain the desired behavior as shown in Figl3] 



B. Example b) n — 3. 

The example we want to analyze corresponds to the case n — 3. Of course we are not going to discuss all the 
possible cases. However n = 3, represents a fair compromise of generalization. In the region where relation (|4H) is 
valid, the integration of Eg. ([^5]) gives 



SnGEj, 167r2aV2 
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FIG. 3: Plot of A/8ttG as a function of the scale invariant x for a = 1/4. The plot works well for backgrounds of the dS, AdS 
and Minkowski type. Note that to obtain a positive induced cosmological constant vanishing at small and large x, we need 
only one parametrization. 



'v^ (2a (-3 + 2x^a) (3 + (-15 - 4,x^a (-3 + x^a)) 7) erf {x^/a) + 



+2e=" °'x^ \-Qal3 - I57 + 2x^07 + 2xaKi j \ -2{a + 25) cosh j + x^af^sinh J^— 

where [x) [v = 0, 1) are the Bessel functions, 5 = P2 and 7 = /33. Instead, in the region where 

mg (r) = 7712 {r) — m\ (r) , Vr € (ri, cx)) , 

we get 

A e-^" 



87rG£;4 87r2a7/2 



{-v^ (157 + 4a;^a2 (/3 + 2:^7) + 6a (/3 + 2a;27)) 



(77) 



(78) 



+2e"^x2a3/2 - x2a<5Xo ( ^) - (4(5 + a (2 + x^^)) i^i 



The asymptotic expansion of Ea. (|77l) in the small x regime is: 



A _ 8a3/2 + 60Fa^ + 157^7 



87r2a7/2 

whereas the leading contributions to Eq. ([77)1 for large x are: 

A _ X (2v^a3/2 + 4a^ + 87 + 3y/^y/aS) G^a^/"^ + Ida^ + 487 + 



327r2a;a4 



(79) 



(80) 
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3 (40v^a3/2 + 128a;9 + 5127 + lOSy^VaJ) 



10247r2a;3a5 



-^0{x-^). 



(81) 



Again, as in the case of n = 1, we find that there is in principle a leading linear divergency in the large x regime. 
However we can choose the parameters satisfying the Minkowski hmit (i.e. the limit of vanishing cosmological constant 
density). This time, differently from the n = 1 case, we can ask that the Minkowski limit is satisfied both in the 
a; — >■ and in the x -> oo region, with a unique choice of the parameters. From Eg. ([50)1 and Eq. ([5T|) follows that the 
parameters that satisfy these requests have to solve the system 



2V7Fa3/2 4q,^ + 87 + 'i^^5 = 
8a3/2 + Q^ajS + 150F7 + = 



(82) 



Notice that already in the case of three parameters (i.e. 7 = 0) the system (15^ can be solved but one gets negative 
values of the cosmological constant density in a large a; > 1 zone (see e.g. FigH]). 




FIG. 4: Plot of A/SttG as a function of the scale invariant x for a = 1/4, /3 = 2y'7r/(97r 
The plot shows that the Minkowski limit is satisfied both in the a; — > and in the x ~ 
present. 



- 32), (5 = (8 - 37r)/(187r - 64), 7 = 0. 
> 00 limit, but a large A < zone is 



Finally, in the full four parameter case (i.e. 7 7^ 0), the system ([5^ can be solved with the further request 
of approaching the Minkowski limit maintaining positive values of A. A solution of Eq. (j82p satisfying this further 
request is: 



a = 1/4, 

136357r + 2048v^ - 38784 



(83) 



1024(97r - 32) 
7687r + 9O90F - 2048 
512(97r-32) ' 



303 
^ ~ ~2048' 

The resulting A as a function of x is plotted is FiglS] A small zone in which the cosmological constant density 
maintains a negative value is however still present. 



IV. SUMMARY AND DISCUSSION 



Motivated by the promising results obtained in the application of Gravity's Rainbow to Black Hole Entropy 
computation [3II and, on the other side in NCG application of ZPE evaluation [l3|. in this paper we have consid- 
ered how Gravity's Rainbow influences the UV behavior of ZPE. We have found that, due to the arbitrariness of 
gi (E/Ep) and 52 (E/Ep), it is always possible to find a form of the rainbow functions in such a way the expression in 
([7]) be UV finite. As introduced in Ref.|5i], the finite result A/SttG is interpreted as an induced cosmological constant 
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A(A-) 




FIG. 5: Plot of A/SttG as a function of the scale invariant x for values of the parameters given by Ea. (|83p . The plot shows 
that the Minkowski limit is satisfied both in the a; — > and in the a; — > oo limit. Both limits are approached maintaining a 
positive A. A small A < zone is however still present. 



but without a regularization and a renormalization to keep under control the UV divergences. To fix ideas we have 
used Gaussian regulators. In this way our approach is directly comparable to NCG . The first evident diflFerence is 
that in NCG the regulator comes into play into the counting of nodes, while in Gravity's Rainbow appear in both 
the sum over eigenvalues and in the counting of nodes. If one fixes the attention on the pure Gaussian regulator one 
discovers that the ZPE is always negative for Gravity's Rainbow. This unpleasant feature can be corrected with the 
introduction of a polynomial with real arbitrary coefficients like in Eq. ()57|) . By imposing the positivity of the result 
for every x = \J (r) / E\, we find that, in case condition (|^ be satisfied, the parametrization ([S7)) is sufficient 
to guarantee that A/SttG > as shown in FigHl On the other hand, when condition (PT|) is satisfied we need two 
different parametrizations to guarantee a correct behavior of A/SttG for x € (0, +oo) and most importantly a point 
of connection where a smooth transition can happen as shown in parametrization ([72]) and in FiglS] In summary, 
the final plot becomesNote that the transition clearly highlights that we need two metrics with the same background 
h (r) but with two different choices of g\ {E/Ep) and g2 (E/Ep). Of course this transition happens because we insist 
to have positivity and a vanishing behavior at the boundary of the range (0, +oo). The vanishing behavior for small 
a; is a guarantee that the Minkowski limit is reproduced. Note that the reproduction of a Minkowski limit in NCG 
(i.e. the — > limit) is less trivial because of the existence of the IR/UV mixing 14j. It also appears that the case 
n — 1 seems to be special because it is the only one that has a positive A/SttG for x G (0, +oo) when condition (|1T|) 
is satisfied. We can focus our attention on this case and suppose to consider the de Sitter background which is the 
static representation of the Friedmann- Robertson- Walker model. For this choice, the shape function b (r) is 

b{r)^^r' (84) 

and the effective masses ([55)1 become on the cosmological throat Tc — Ads 

"^0 (r) = ml (r) = (r) = Ads- (85) 

Since the behavior of A/SttG for the de Sitter universe is described by Fig (31 we can compare our results with 
observation. Since A/SttG represents the observed cosmological constant induced by quantum fluctuations of the pure 
gravitational field, we can fix its value at the present day as 

8^ - 10-"ey^ (86) 

which can be described either for x ^ 1 or a; ^ 1. However, for the de Sitter case x ~ \fAds/Ep and x means 
Tc ^ 1 which is in agreement with our present universe. On the other hand, when a: 3> 1 means that Tc ^ 1 which 
could be in agreement with the very early universe except for the disagreement with the expected theoretical prediction 
which for our plot in units of Ep should be O (1). Therefore it appears that only the left branch of Fig 13] from the 
bottom to the hilltop can be interpreted as a sort of a backward evolution" in the radial coordinate r. However to 



14 




FIG. 6: Plot of A/SttG as a function of the scale invariant x for values of the parameters given by Ea. (|72p . The plot shows that 
the Minkowski limit is satisfied at the boundaries of the range (0, +oo). Both limits are approached maintaining a positive A. 
The undesired part of the plot has been eliminated to visualize the global behavior. 



follow the curve from x ~ 2.18 to a; ~ one should have a variable A^s, namely Kds = ^ds {f)- The same situation 
appears to exist for Minkowski space in radial coordinates and for the AdS space. Fortunately, Minkowski space has 
no a preferred scale and Figl3] has the correct asymptotic behavior except for an unpleasant peak in correspondence 
of the peak location of the dS space. It is likely that this spurious prediction is due to the coordinate choice. On the 
other side one can verify that A/8ttG — when r — > for both Minkowski, dS and AdS spaces. Coming back on the 
AdS space, we have to note that this background is not endowed of a horizon and therefore it makes difficult to find a 
significant point. Nevertheless looking once again Figl21 we can claim that for r — >■ oo and very small AAdS one gets 



6(r) 



^AdS 3 

:z — r 



and 



l6/r^ + AAdS 
El 



^0, 



(87) 



namely a vanishing K/S'kG, while in the other regime, i.e. r — > 0, a; — >■ oo and once again one obtains a vanishing 
A/BttG that it means that against all odds, we have regularity on the singularity r = 0. 

Appendix A: Integrals 

In this appendix, we explicitly compute the integrals coming from Ea. (|^5)) . We begin with 



/+ = 3 



" E^ 
^C— exp(-a-^; 



E^dE'^+ml {r)dE. 



(Al) 



It is useful to divide /+ into two pieces with i odd and i even, thus we can write 



(A2) 
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where 



4 = Z/2E% c. (-)' lim — / exp [-{a + Jx + (r) dx, 

II = ^/2E% V Ci {-r lim — / a; exp [-{a + I3)x] Jx + ml (r) /Eldx, 
i=o •'O 

are expressed in terms of the variable x = /Ep. 

The integrals involved in the expressions of and /° can be evaluated using the formulas 

^ dx{x + t)^^^x^^'^exp{-i^x) = ^exp(^^^ Ki(^^^ i>0,/i>0 

/ dx{x + t)^^^xexp{-i^x) = -^ + ^IJ-^/^exp{tii){3-2ti^)Eik[y/t^i\ t>0,fi>Q 
Jo 2^4 



obtaining 



/-3/2i.^^c.(-riim- 



i=0 



ml (r) 



exp 



4 (r) (a + /3) 



ml (r) 

2£;|, 



/; = +3/2i^^^c,(-riim- 



■ exp 



i=0 
2 



3mo (r) 



TO^ (r) (a + /3) 



2ml ('') 



^|,(a + ^)3/2 (a + /3)5/2 



Erfc 



— — v(" + P) 



The same procedure can be followed to evaluate 



^a-^expi-a-^) 
Li=o P 



E'^JE^-ml {r)dE. 



Even in this case, it is useful to divide /_ into two pieces with i odd and i even. Thus, we can write 

/- = 11+11, 

where /f. and 7° are given by 



It = 3/2Ej, J2 (-)' Toi / + - ™o (0 /E^dx, 

^-^/^^^E^^HMiml- 
Using now the formulas 



i=0 



ap« Jm^/El 



.1/2 



exp [—(a + P)x] '^x — ml{r) / Epdx. 



dx{x-tf''^x^l'^exp{-iJLx) = ^exp(^-^-^Kx(^-^ t>Q,ii>0 
dx{x - t)^^"^ xexp{-ij,x) = ^n~^/^{3 + 2ijt)exp{-ij,t) t > 0,/i > 



(A3) 
(A4) 



(A5) 
(A6) 

(A7) 



(A8) 

(A9) 

(AlO) 

(All) 
(A12) 

(A13) 
(A14) 



It and 7° can be rewritten in the form 



3/2Eiy"ci(-y lim 

^ /3-s-O 



ml (r) 



dpi [2{a + P)Ej, 
Il=3/2Eip,i-)Jrn^^^[^ia + py 



ml (r) (a + /?) 



27;|, 



Ki 



ml (r) (a + /3) 



27;|, 



5/2 



3 + 2(a + /3)4W 



exp 



ml (r) (a + /3) 



El 



(A15) 
(A16) 
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